clc;
clear;
f=@(x)4*x(1)*x(3)*exp(-x(1)^2*x(2)-x(3)^2);
tic,
I=quadndg(f,[0 0 0],[2,pi,pi]),
toc